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ABSTRACT 



The average mass density profile measured in the CNOC cluster survey is well 
described with the analytic form p(r) = Ar~ l {r + Ob p ) , ctS advocated on the basis 
on n-body simulations by Navarro, Frenk & White. The predicted core radii are 
a p = 0.20 (in units of the radius where the mean interior density is 200 times the 
critical density) for an £1 = 0.2 open CDM model, or a p = 0.26 for a flat O, = 0.2 model, 
with little dependence on other cosmological parameters for simulations normalized 
to the observed cluster abundance. The dynamically derived local mass-to-light ratio, 
which has little radial variation, converts the observed light profile to a mass profile. 
We find that the scale radius of the mass distribution, 0.20 < a p < 0.30 (depending on 
modeling details, with a 95% confidence range of 0.12 — 0.50), is completely consistent 
with the predicted values. Moreover, the profiles and total masses of the clusters as 
individuals can be acceptably predicted from the cluster RMS line-of-sight velocity 
dispersion alone. This is strong support of the hierarchical clustering theory for the 
formation of galaxy clusters in a cool, collisionless, dark matter dominated universe. 
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1. Introduction 

A fundamental prediction of the hierarchical clustering paradigm in a collisionless dark 
matter dominated universe is that the average mass profile of dark halos must primarily depend 
on mass, if the density perturbation spectrum is nearly a power law. The single parameter surface 
brightness profile, £ oc exp (— r 1 / 4 ), was proposed by deVaucouleurs (1948) to describe elliptical 
galaxies, but it was later found also to be quite an accurate description of isolated n-body collapse 
simulations, both with and without dissipation ( van Albada 1982 , |Carlberg, Lake fc Norman 



1986). This r 1//4 "law" has no simple analytic form to describe the volume distribution ( |Young| 



1976), motivating Hernquist (1990) to propose the function ph{t) = Ar l (r + a) 3 as a more 



tractable alternative, which in projection is usefully close to the r 1//4 function. 

In cosmological halo formation simulations, where external tides and continuing infall are 
included, the central density profile is found to remain approximately r _1 (Dubinski &; Carlberg 



1991, Carlberg 1994, [5yer fc White 1996j ), although there is growing evidence that the central 



slope is somewhat steeper (Fukushige & Makino 1997, Evans & Collett 1997). However, beginning 



near the virialization radius the density profile is found to fall somewhat less steeply than in 
isolated collapses, leading Navarro, Frenk & White (1996 and 1997, hereafter NFW) to advocate 
p(r) = Ar _1 (r + a)~ 2 as an empirical description of the halo density profiles. In a series of 
exhaustive n-body simulations the relation between the scale radius, a, of the halo and its mass 
has been studied to yield an empirical understanding of its value ( |Cole Sz Lacey 1996 , NFWj ). 



The principal result of interest here is that the scale radius, a, is measured in the simulations to 
be about 20% of the "virial radius" for the 10 15 M Q dark halos appropriate to rich clusters. This 
result has only a weak dependence on either the power spectrum, or the world model parameters 
Q and A, for simulations normalized to the observed cluster abundance. Halos of much lower 
mass, say the 10 12 Mq characteristic of individual galaxies, are generally much more centrally 
concentrated than clusters, with quite a strong dependence on the cosmological details. 

In the next section we fit the galaxy number density profile, v{r), to the NFW function, then 
present evidence that this profile (specifically for galaxies selected in Gunn r, which as a red pass 
band is relatively insensitive to the current star formation rate) accurately traces the average 
cluster mass profile, p{r). In particular we compare the derived surface mass profile with the 
surface galaxy number profile. Our analysis is independent of the light-traces-mass assumption 
and allows for a range of velocity dispersion profiles. We have previously shown that this type of 
analysis works to recover the same mass profile from two vastly different subsamples of the data 
(Carlberg et al. 1997a). The fitted scale radius of the NFW model is compared to the predicted 



scale radius. In Section 3, we use the scaling from velocity dispersion to radius to predict the 
form of the profile for each cluster in turn and compare it to the available data. The results are 
briefly discussed in the final section. We use H = 100 km s" 1 Mpc" 1 and ft = 2<? = 0.2 for all 
our calculations, although the results are not very sensitive to these choices. 
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2. The Average Mass Profile of a Galaxy Cluster 



The Canadian Network for Observational Cosmology (CNOC) Cluster Survey (Carlberg 
et al. 1994, Yee, Ellingson fc Carlberg 1996j ) obtained ~ 2600 redshifts of Gunn r selected galaxies 
in the fields of 16 high luminosity X-ray clusters at z ~ 73. This is a relatively homogeneous 
sample of clusters that are guaranteed to be at least partially virialized on the basis of their X-ray 
emission. The virial mass-to-light ratios of these clusters were found to be identical, within their 
measurement errors of ~ 25% ( Carlberg et al. 1996 ). The 14 clusters that are "non-binary" are 



combined in normalized co-ordinates to make an ensemble cluster. This diminishes substructure 
and asphericity to a level where the galaxies can be treated as if in a spherical distribution which 
is consistent with dynamical equilibrium ( Parlberg et al. 1997a| ). Galaxies of all colors above a 
k-corrected Gunn r absolute magnitude of —18.5 are included in the average. 



2.1. The Average Cluster Number Density Profile 

To combine the clusters the brightest cluster galaxy (BCG) is used as the nominal center of 
each cluster on the sky following our earlier procedures (Carlberg, Yee & Ellingson 1997, hereafter 



CYE). The galaxy velocities are normalized to a\, the RMS velocity dispersion of the clusters 
about the cluster mean. The projected radii are normalized to an empirically determined r^oo; the 
radius where the mean interior overdensity is 200/9 c . To derive r2oo from the observational virial 
radius, r v (which is largely fixed by the outer boundary of the sample), we assume that M(r) oc r. 
This gives 

7-200 - To W)' { ) 

which is completely independent of the observational virial radius. Usually the extrapolation in 
radius is a modest 25%, so the precise mass profile assumed makes little difference to the result 



flCYED 



The average projected galaxy number density profile, S(i?), is fit with the projection of the 
volume density function, 

^ = ( t (2) 
r(r + a u ) p 

where p is fixed at either p = 2 (NFW) or p = 3 (Hernquist). The results for p = 2 are shown in 
Figure || as the solid line. The dashed lines are described below. Both the p = 2 and p = 3 forms 
are statistically acceptable fits by the x 2 test. We will only consider the p = 2 form for the rest of 
this paper. The fitted scale radius is a u = 0.27 with a 95% confidence range of 0.13 — 0.43 from 
the x 2 distribution. 
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2.2. The Relationship between the Mass and Number Density Profiles 

There is no dynamical necessity for v(r) to be directly proportional to p{r). The relationship 
between the two is derived from the projected velocity dispersion, a p (R). The dynamical mass 
profile, Mjy(r), which is the volume integral of p(r), is inferred from the Jeans Equation (e.g. , 
Binney fc Tremaine 19871 , |CYEj ), 



air 



M D {r) = —fr 



d In al d\nv 
a In r amr 



(3) 



where the velocity anisotropy parameter, (3 = 1 — oj/cr 2 . N-body simulations for a variety of 
cosmologies show that the dependence of (3 has a nearly universal radial variation (|Cole fc Lacey 



1996) which we somewhat arbitrarily model as, 



P(r)=P m ^—, (4) 

which takes on a maximum value of (3 m at r = 2, in r2oo units. We demonstrate below that the 
volume mass density is not very sensitive to the details of the velocity modeling, provided the 
model is an acceptable statistical description of the data. We use values of 0.3 and 0.5 for the 
parameter /3 m which roughly bracket the range of (5 seen in the simulations. 

We model the radial velocity dispersion as, 

2 _ D cir/(1 + Ci r) +c 2 

where B and b are the two parameters adjusted to fit the observed cr p (R). We include here a new 
central data point at 0.05r2oo- The rest of the data are as in CYE| . The [ci,C2] parameters are 



externally fixed to allow us to vary the shape of the curve. In Figure [2] we display results where 
[ci, C2] are [0, 1] and [1, 0] to give the two extremes of the functional form. The [0, 1] form was used 
in |CYE| and [1,0] gives a dispersion which drops to zero at the center. The intermediate values 
[8, 1/2] are adopted for reasons made clear below. Both the [1, 0] and [8, 1/2] forms have reduced 
X 2 near unity. The [0, 1] form is about 4% probable, although this is critically dependent on the 
0.05r2oo point, the concern being that the old central galaxies may be a distinct, low velocity 
dispersion, cluster population. 

We show the results of the mass analysis in Figure |3[ The quantity plotted is p(r)/v(r) 
normalized to M v /L x L(r), computed from the same data. The ratio p(r)/u(r) is only weakly 
constrained inside 0.1r2oo an d the assumption of virialization likely fails outside 1.5r2oo- Within 
this radial range p(r)/u(r) varies less than about 30% of its mean value. It should be noted that 
for most purposes the more relevant quantity is Mp^r), the volume integral of p(r), which is 
substantially more slowly varying. The offset of p(r)/v(r) below unity indicates that mass and 
light are similarly distributed but that the virial mass has a scale error. A wide range of positive 
intermediate values of c\ and ci gives nearly constant buv everywhere. We display the c\ = 8, 
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C2 = 1/2 case, which is one of many intermediate sets of parameters which leads to substantial 
cancellation of the opposite behaviour of the two extreme a r (r) models in the derived Mo{r). 

With the dynamically derived p{r)/v(r) in hand we can use it to model the galaxy number 
profile, ^(r), assuming the mass profile is described by the NFW profile. This procedure measures 
o p , the scale radius of the mass distribution. We take p(r)/u(r) to be equal to the endpoint of 
the range for r < 0.1r2oo and r > 1.5r2oo- The results for our three velocity dispersion functions 
and the two values of /3 m are shown in Figure [l] as the dashed lines. The resulting a p values range 
from 0.20 ([ci,C2] = [0,1]) to 0.30 ([1,0]). This expands the 95% confidence range of a p over a u 
slightly, to 0.12 < a p < 0.50. We conclude that galaxies selected from their red light brightness, 
i.e. , a passband well to the red of the 4000A break to minimize sensitivity to star formation (and 
not necessarily galaxies red in color, parlberg et al. 1997a| ), accurately trace the mass profile 
of the cluster. That is, u{r) is directly proportional to p(r), within the statistical limits of our 
measurements. 

3. Individual Cluster Profiles 

The NFW profile can be tested for its application to the clusters as individuals. The r2oo f° r 
each cluster is derived from the observed velocity dispersion, Equation [l]. The NFW scale radius 
is set at our best fit value, a u = 0.27. The fits are restricted to the region inside 1.5r2oo> which 
is expected to be virialized for Q ~ 0.2. The resulting x 2 P er degree of freedom (there being 
10-16 degrees of freedom) are given in Table [lj. The cluster MS0906+11, which we observed, is 
excluded because it was not possible to measure a reliable velocity dispersion, which we believe 
is likely because there is another nearby clusters in the redshift direction (|Carlberg et al. 1996| ). 
The reduced x 2 values of 14 of the 15 clusters are 1.1 or less, indicating that the NFW function 
with a fixed scale radius is a statistically acceptable fit to the density profile. The A2390 cluster 
has x 2 = 1-38 per degree of freedom which indicates that the fit is about 15% probable. A2390 is 
the cluster with the most data. If we had about a third less data for this cluster, the fit would 
have been entirely acceptable. However, as the amount of cluster data increases beyond about 
100 galaxies, the substructure within the cluster ( [Abraham et al. 1996| ) becomes more clearly 
detected. By extension, if we had obtained more than 200 or so galaxies within any cluster, the 
substructure would have become clearly resolved. This sampling strategy was part of the initial 
layout of the CNOC observations. 

4. Discussion and Conclusions 

We have established that the dynamical mass profile, p(r), is proportional to the galaxy 
number density profile, i^(r), over the range of radii where the argument is secure. We further 
make the modest assumption that v{r) continues to trace the mass, within the errors, at larger 
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radii. Consequently we conclude that the NFW function accurately describes the mass field. 
The measured scale radius is 0.20 < a p < 0.30 with a 95% confidence range of 0.12 — 0.50. 
This relatively large confidence range is a result of the relatively small change in profile slope 
at the characteristic scale for the NFW profile. The Hernquist function has much smaller errors 
when used in same fitting procedure. The predicted values for Q = 0.2 and a CDM spectral 
normalization of as = 0.95 are 0.20 for an open model (nearly identical to the as = 0.6, 0, = 1 
value of 0.19) and 0.26 for a flat low density model ( NFW|) , are remarkably precise predictions of 
the observed value. 

Overall it appears that a very good understanding of the evolution of the mass field in 



clusters is emerging. The value of 0, determined from dynamical observations of clusters ( Car lb erg 



et al. 1996, CYE) predicts a slow change in the abundance of clusters with redshift, which is 



in accord with the directly observed value ( |Carlberg et al. 1997b ). Previously we found that 



clusters as individuals have masses inside r2oo which are accurately predicted from their velocity 
dispersions ( Parlberg et al. 1996 ). Here we have shown that the predicted scale radius and 



the universal NFW profile provide a statistically acceptable description of the virialized region, 
r < 1.5r2oo 5 of 14 of the 16 clusters selected for the CNOC observational program, with the 
failures likely being the result of large substructures. The combination of successful predictions of 
galaxy cluster mass profiles and their cosmological number density evolution is a strong argument 
for the model of cluster formation in an effectively cold, collisionless, dark matter dominated 
universe. These results could be substantially tightened, in the first instance, through redshift 
determinations for more cluster galaxies in the outskirts of this sample. 



The referee, Cedric Lacey, provided valuable comments which substantially improved the 
presentation of this paper. We thank CFHT for the technical support which made these 
observations feasible. NSERC and the NRC provided financial support. 
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Table 1: Goodness of Fit of A Universal Profile 
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Fig. 1. — The average galaxy density profile of the clusters using the BCG as the center. In the 
text, we demonstrate that this galaxy number density profile is statistically equal to the projected 
average mass profile. The solid line is the projection of the NFW function. The dotted lines are 
the projection of the product of an assumed NFW mass profile with the dynamically derived local 
light-to-mass ratio, u(r)/p(r). The dashed lines are the projected NFW mass profiles themselves, 
offset downward by half a dex. 

Fig. 2. — The projected velocity dispersion profile and the projection of the fitted profiles, for 
both (3 m = 0.3 and 0.5 which are nearly indistinguishable on this plot. The dotted line is for 
[ci, C2] = [0, 1], the dashed for [1, 0] and the solid is [8, 1/2] (see text for details). 

Fig. 3. — The derived ratio of the dynamical mass density profile, p(r), to r selected galaxy profile, 
v(r), normalized with the virial mass-to- light ratio evaluated inside 1.5r2oo- I n each pair of curves 
the upper line at small radius is for (3 m = 0.3 and the lower for (3 m = 0.5. The dotted line is for 
c\ = 0, the dashed for C2 = and the solid is our preferred C\ = 8, C2 = 1/2. 




Fig. 1.- 



- 11 - 




- 12 - 




0.5 1 1.5 2 

r / r 200 



Fig. 3.- 



